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Large nuclear Overhauser fields detected in vertically-coupled double quantum dots 
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We report the electrical induction and detection of dynamic nuclear polarization in the spin- 
blockade regime of double GaAs vertical quantum dots. The nuclear Overhauser field is measured 
within a certain range, relying on bias voltage control of the interdot spin exchange coupling and 
measurement of dc current at variable external magnetic fields. The largest Overhauser field ob- 
served was about 4 Tesla, corresponding to a nuclear polarization ~ 40% for the electronic g-factor 
typical of these devices, \g*\ w 0.25. The amount of polarization that can be detected using the 
present methods is limited by the range of externally applied magnetic fields over which a clear cur- 
rent step can be observed. A phenomenological model based on the hyperfine mixing of electronic 
singlet and triplet spin states is proposed to explain these observations. 

PACS numbers: 73.63.Kv, 72.25.Rb, 76.70.Fz 



Abundant nuclear spins present in semiconductor 
nanostructures couple collectively to the spins of con- 
fined electrons through the contact hyperfine interaction. 
In situations of mesoscopic confinement, an electron spin 
typically interacts with TV ~ 10 5 — 10 6 nuclear spins. Gal- 
lium arsenide (GaAs) consists of 71 Ga, 69 Ga and 75 As iso- 
topes, each bearing nuclear spin I = 3/2. The hyperfine 
coupling constant in GaAs is A sw 90/ieV [lj. For com- 
pletely polarized nuclear spins, this interaction leads to 
nuclear Overhauser fields ~5-10 Tesla, depending on the 
precise value of the electron g-factor in the nanostructure 
[l|. Such large nuclear fields can produce measurable ef- 
fects on spin-dependent electron transport; for example, 
in the spin-blockade (SB) regime of double quantum dots 
@> El or in quantum well structures in the quantum Hall 
regime at special Landau level filling factors [4j. In the 
absence of near-perfect nuclear polarization, field fluctu- 
ations order A/yN are the main source of dephasing for 
electron spin qubits in GaAs quantum dots B], 0, 0, [I| • 
Dynamically polarizing the nuclear spin system either op- 
tically 0, [l(| or electrically Elhas been proposed as one 



way to mitigate this effect [5j, [12[ . It was also proposed 
that polarized ensemble of nuclear spins could serve as 
qubits or quantum memory, taking advantage of radio- 
frequency manipulation, transfer of information to elec- 
trons, and very long coherence times 1^, 14 1. Developing 
electrical control of the nuclear polarization in such de- 
vices would therefore be an important contribution to the 
present suite of experimental controls. 

Earlier experiments on transport in the SB regime of 
vertical coupled dots showed magnetic-field dependent 
current features such as a hysteretic current step and 
slow oscillations Q. These features were identified with 
the presence of dynamic nuclear polarization (DNP), as 
well as complex feedback between the electronic and nu- 



clear spin systems. It has been proposed that DNP oc- 
curs in this system near a level degeneracy at which a 
blockaded spin-triplet state with one electron in each 
dot (e.g. | XI) = | J.J.)) is mixed with the singlet state 
\S) = 4=(| U) - I IT)) via the flip-flop terms of the hy- 
perfine interaction, allowing transport through the dots 
[2, EH. A mutual electron/nuclear spin flip/flop occurs 
to conserve angular momentum, and energy conservation 
is in general provided by inelastic processes Since 
the hyperfine mixing is only efficient for one of the triplet 
states, the nuclear spins are pumped toward a polarized 
state. Cotunneling processes due to strong coupling with 
the leads yield finite lifetimes for the triplet states not 
efficiently mixed by the hyperfine interactions. The accu- 
mulated nuclear polarization creates an Overhauser field 
(henceforth defined as the average Overhauser field of the 
two dots) that shifts the energy levels of the \T±) triplet 
states. In this Letter, we determine the Overhauser field 
based on dc current measurement and voltage control of 
the interdot exchange coupling, with results that imply 
large nuclear polarizations. The degree of nuclear polar- 
ization can be controlled by varying the external mag- 
netic field. Finally, a phenomenological model is pro- 
posed to provide a qualitative understanding of these re- 
sults. 

Experiments were performed on a_0.4//m diameter ver- 

HEl 



tical double dot structure 



lal with 10 nm GaAs 



quantum wells, 7 nm Alo.3Gao.7As outer tunnel barriers 
and 6.5 nm center barrier, at a sample temperature 1.7 
K. Gate voltage was adjusted to give a clear SB region 
in the current-voltage characteristic 0, EH- Measure- 
ments consisted of recording dc current as a function of 
source-drain voltage at variable external magnetic fields. 
Crucially, the external field was applied in the plane of 
the two-dimensional electron gas (i.e. the lateral plane), 
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FIG. 1: (upper, left) Scanning electron microscope image of a 
device nominally identical to the one investigated here. The 
double quantum dot is located in the pillar structure just be- 
low the surface of the gate metal, (upper, right) Spin-blockade 
region (labeled 'SB') in the I-V characteristic at zero external 
magnetic field (B ex t = 0). (lower) Schematic energy diagram 
of the two-electron eigenstate energies as a function of the rel- 
ative detuning e of the two dots. The main figure shows the 
anti-crossing (resonance) of the (1,1) and (0,2) charge config- 
urations of the spin singlet state, which are coupled by the 
interdot tunneling parameter t. The (1,1) triplet state energy 
is independent of detuning. The expanded view shows the 
region just to the right of the resonance which is the focus 
of these experiments. The Zeeman splitting of the T+/T- 
states depends both on the external field B ext and the aver- 
age nuclear Overhauser field B n . The spin exchange energy 
is defined as J, and the level separation between |T_) and \S) 
is denoted by E S /t_ ■ 



so that it had negligible effect on the dot electronic wave- 
functions. 

Figure [1] shows a schematic energy level diagram 
for the two-electron states relevant to these experiments. 
The horizontal axis is the degree of energy detuning e 
between the two dots, which is varied experimentally by 
changing the source-drain voltage V s d- In our device, 
e » Q.TJVsd + e o, with e = located slightly to posi- 
tive V s d, however the precise position is difficult to de- 
termine from analysis of the standard Coulomb diamond 
measurement. The energy eigenvalues plotted in figure Q] 
were calculated by taking into account intra- and inter- 
dot Coulomb energies, the interdot tunnel coupling t, and 
e(V s d) 15]. An interdot tunneling resonance between the 
(1,1) and (0,2) charge configurations of the spin singlet 
state is centered at e — 0. At large negative detuning, 
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FIG. 2: (a) Expanded view of the spin blockade region, show- 
ing current steps obtained by sweeping V s d (left to right) at 
B ext = 0.2 T and B ex t = 2.4 T. V s te P indicates the position 
of the 0.2 T step. A dot-lead resonance peak near V a d = 1 
mV (see upper right of figure 1) gives rise to the background 
slope for Vsd < 2.0 mV. (b) 'Relaxed' trace is the B ext = 0.2 
T trace in o; 'Polarized' trace was obtained at the same exter- 
nal field after pumping the nuclear polarization at Vsd = 3.0 
mV for 30s. 



the (0,2) singlet is high in energy and Coulomb block- 
ade occurs. Spin-blockade is mainly observed at posi- 
tive detuning (i.e. to the right of the resonance). The 
spin exchange coupling J is defined as the energy differ- 
ence between the singlet and spin-zero triplet |Tn), with 
a maximum value J| £= o = s/2t. The expanded region 
in figure [T] details the region of interest for our experi- 
ments. It shows that application of an external magnetic 
field B ext splits off the \T±) states due to the Zeeman 
interaction. Nuclear polarization arising in the two dots 
from the hyperfinc mixing of |T_) — > \S) is antiparallcl to 
B ext , but causes an average Overhauser field B n || B ext 
since sgn(gA) = —1. With respect to the \T±) energies, 
B n and B ext are equivalent; this fact is central to our 
method for detecting B n . Note that the energy differ- 
ence E S / T _ between |T_) and \S) plays a leading role in 
the DNP process under study. By varying detuning us- 
ing V s d, we can manipulate the field B ext + B n for which 
E s/T _ = 0. 

Figure [2k shows that a current step of magnitude 
~0.5 pA is observed when V s d is swept towards larger 
detuning at fixed external field B ext [15j | . The character- 
istic step is very similar to that observed in Ref. Q while 
sweeping B ext at fixed V s d- This step appears to indi- 
cate the level crossing E$/T- — 0; an d it shifts toward 
smaller V s d (i.e. smaller detuning) with increasing B ex t 
as expected from figure [TJ We assign the step position 
Vstep to the edge of the upper current level, as shown in 
figure^. No step was detected for B ext < 0.1T in the 
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FIG. 3: (a) Measured current step positions Vstep versus ex- 
ternal field, for the two measurement sequences described in 
the text. The data from the 'relaxed' sequence have been nu- 
merically fit to an empirical function. The differences along 
the vertical axis between the 'polarized' data points and this 
curve yield estimates of the Overhauser fields B„, which are 
plotted in b. (b) The values that rely on the extrapolated por- 
tion of the fitting curve in a (e.g. for B ext > 3.6 T) are shown 
as open circles, and those that can be determined directly 
from the data points are shown as filled squares. Experimen- 
tal uncertainty is indicated by the size of the symbols. 



present device. Traces were obtained at a sweep rate of 
0.15 mVs- 1 . 

Figure [2Jd shows that V ste p is also shifted in the pres- 
ence of a nuclear Overhauser field. To obtain the trace 
labeled 'polarized', we first pumped the nuclear polar- 
ization by fixing V s d > V s t ep for the V s t ep observed in 
the absence of an Overhauser field (trace labeled 're- 
laxed'). After pumping about 30s, V s d is rapidly returned 
to Vo = 1 mV and subsequently swept up again to obtain 
the second trace. Both traces in figure^ are taken at the 
same value of B ext = 0.2 T. The position of the shifted 
current step is nearly identical to the V s t ep obtained at 
larger B ext = 2.4 T in figure [2^. From these experi- 
ments we conclude that the Overhauser field B n present 
in figure [2jj can be estimated as the difference in external 
fields applied in figure [2^. All such measurements were 
preceded by suitably long wait times at V s d — to allow 
for relaxation of residual nuclear polarization from previ- 
ous traces. Relaxation of nuclear polarization occurring 
during the readout sweep is discussed below. 

The experiments of figure[2]were repeated over a range 
of B ext values from 0.1 — 3.2 T. The results are summa- 
rized in figure [3J The points in figure [3^, are the mea- 
sured values of V s t ep versus B ext for the two measurement 
sequences described above. The estimated Overhauser 
fields B n correspond to the differences between the two 
curves along the B ext axis. The B n estimates plotted in 
figure [3Jd were obtained by subtracting the points from 
the 'polarized' sequence from an empirical curve fit to the 
'relaxed' data points. Note that figure [3Jd differentiates 
between B n values that can be determined directly from 
measured data points and those determined from the ex- 
trapolated portion of the fitted curve at B ext > 3.6 T 
in figure [3^. The behavior observed in this latter range 




FIG. 4: A schematic diagram showing the dynamics of the 
|T_) energy level as the source-drain voltage V s d is swept from 
left to right and dynamic nuclear polarization occurs near the 
E,s/t_ — degeneracy. J step indicates the value of J at the 
crossing point, which is approximately the position (in V s d) 
of the observed current step. 



is most likely due to the polarization dynamics occur- 
ring near the zero detuning point, i.e. where the states 
|T_) and \T + ) can both be nearly degenerate with the 
upper and lower singlet branches, respectively. In fig- 
ure[3p, B n rises rapidly in the range 0.1T< B ext < 0.4T 
reaching its maximum value ~ 4 T near B ext = 0.65 T. 
It then falls off quasi-linearly to ~ at B ext = 3.2 T 
(notably, the sum B n + B ext is roughly constant in this 
range 0.65T< B ext < 3.2T). 

Relaxation measurements of the quantity B n (r) were 
obtained by inserting a time delay r at V s d = after 
pumping in the 'polarized' measurement sequence. We 
observe a bi-exponential decay curve with approximate 
decay times of 13 ± 3 and 68 ± 10 seconds [2fjj. In gen- 
eral, this decay will be due to both diffusion of nuclear 
polarization away from the dot region and the intrinsic 
nuclear relaxation processes; however we expect intrinsic 
relaxation to be much slower than the rate of loss due to 
spin diffusion, and we estimate the latter to yield a ~10- 
20s timescale (i.e. for diffusion in the vertical direction). 

Given the finite duration of the readout voltage sweep 
in this scheme, the true Overhauser fields should be 
slightly larger than these values, by a factor ~1.2 for the 
largest fields, if a correction is made based on the mea- 
sured relaxation. Using the expression for B n in GaAs 
obtained in Ref. [l|, inserting |<?*| = 0.25 from measure- 



ment of a nominally identical device 21(, and assuming 
equal isotopic polarizations, we obtain an average nu- 
clear polarization P _B„/(9.32T). The data therefore 
suggest a maximum polarization lying within the range 
0.38-0.52. 

To model this behavior, we consider the rate for an 
electron to escape the spin blockade due to hypcrfinc- 
induced singlet-triplet mixing. The two-electron spin 
Hamiltonian Ti = Tio + Ti' consists of the Zeeman and 
exchange terms (Tio) an d perturbing hyperfine coupling 
terms (Ti.'). The triplet and singlet eigenstates of Ho are 
mixed only by the hyperfine terms antisymmetric with 
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respect to the two dots, i.e. the terms involving dif- 
ferences of their effective nuclear fields. Let us define 
the antisymmetric nuclear fields h a — — ((1^) ~ (-0)> 
where (J^) is the average a component of spin per nu- 
cleus in dot f (here we assume uniform hyperfine cou- 
plings for simplicity). Statistical fluctuations give rise 
to fields \h a \ ~ A/yN where N is the average number 
of nuclear spins. The rate of nuclear spin Am z = — 1 
transitions (T~ ►!) is determined by the rate of escape 
from |T_) via hyperfine mixing with |5(1,1)), and sub- 
sequent tunneling to the (0,2) and (0,1) charge config- 
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this can be treated by 



second-order perturbation theory [1 61 ] and the rate writ- 
ten as = p TL | (5(1, 1)| \T'_ }| 2 r s(M) ^ (0il) , where 

Pt> is the steady-state occupation probability for the 
perturbed eigenstate \T'_) and fsn^ 



♦(o,i! 



is the rate 



of escape from |5(1,1)) to (0,1) via (0,2). In gen- 
eral, Tgn .i)_»(o.i) can reflect both elastic and inelastic 
processes, and therefore depends on tunnel couplings, 
phonon scattering rates and cotunneling processes; here 
it is only considered phenomenologically as a parame- 
ter. First-order perturbation theory gives (S(l, 1)\\T'_) = 
h+Mhi A Therefore W* i w a( A ^ ) 2 

V2E S/T _ ^NE S/T _- lnereIOre ^T^l ~ a \E s/T _) > 

where a = Pt' T 's(i,i)->(o,1) (the opposite rate Wx^j is 
computed the same way using the \T' + ) state). 

The rate of change of nuclear polarization in the j th 
dot can be expressed: 



d{p z ) _ wwi -y 3/2 ) - Ww(i-pL a/2 ) 



dt 
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where p±s/2 are the normalized populations of the m z = 
±3/2 sublevels and we assume an equal probability 1/2 
for a nuclear spin transition to occur in dot 1 or dot 
2. Tx is the phenomenological relaxation rate of nuclear 
polarization (including spin diffusion out of the dots) 
and Nj is the number of nuclear spins in the j th dot. 
As V s d is swept up, Egj T _ decreases (see figure Q} and 
Wf_ >| increases; a small nuclear polarization begins to 
accumulate. Since E S / T _ = J — \g*\pBB ext + A{Iz), 
Es/T- is further decreased by the nuclear Overhauser 
field A(I Z ) = —\g*\/J>BB n < and positive feedback 
occurs. Once this occurs, the system is forced to pass 
through the degeneracy E S / T _ — where the pump- 
ing rate is maximum and is limited by the (non-nuclear 
spin- flip) escape rates of the \T + ) and \Tq) states [22j |. 
This is illustrated schematically in figure HI A steady- 
state is then reached in the regime of E s / T _ < 0, at 
which the steady-state nuclear polarization can be deter- 
mined by setting d ^^ ss = and solving self-consistently 
for (Iz)ss- Assuming H^_>x >> and setting 

P±3/2 — 1/4 (low polarization limit), we find for the 
steady-state polarization: 
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to this analytical solution but would be determined by 
a full simulation of the dynamics. We take equation [3] 
to be valid for < [3 < 1/2 Inserting T x = 30s, 

N = 5 x 10 5 and (3 = 1/2, equation [2] requires that 
a - 1.5 x 10 8 Hz for \P SS \ - 20%. Such a fast rate for 
a (on the order of the dot-lead tunneling rate) is possi- 
ble near the singlet-singlet resonance where 15(1, 1)) and 
|5(0, 2)) are hybridized by a relatively large interdot tun- 
nel coupling ~ 50^ieV. A full explanation of the depen- 
dence of B n as a function of B ext , and of the current step 
24} , is beyond the scope of this paper, and will require a 
detailed modeling of the double-dot transport and polar- 
ization dynamics. However, we do expect B n to fall off 
as B ext is increased (as in figure [3j3 for B ext > IT) due 
to increasing pumping of the opposite transition VFj^j 
via the state |T+) mixing with the lower singlet branch. 

In conclusion, we have shown that surprisingly large 
nuclear polarizations can be electrically induced and de- 
tected in the present GaAs double dot system. One draw- 
back is the detection scheme: the magnitude of polariza- 
tion that can be detected is limited by the range of B ext 
over which we can observe a clear current step, and ad- 
ditionally by the fact that dot-lead tunneling rates will 
be reduced at large B ext due to a momentum mismatch 
between dot and lead states [lj| . Indeed, it is likely that 
significantly larger polarization was actually achieved in 
our experiments than was detected, since the current step 
observed in the readout V s d sweep suggests a secondary 
polarization boost. Another detection method, such as 
electron spin resonance, could be used to extend these 
measurements in similar devices. 
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